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ABSTRACT 

A  theory  is  developed  which  describes  the  local  wave 
effects  near  a  wave  front  advancing  into  calm  water  and  is 
uniformly  valid  near  the  front.   The  waves  generated  by  an 
oscillating  vertical  boundary  (a  wavemaker)  and  propagating 
out  into  an  infinitely  long  channel  are  investigated  both 
theoretically  and  experimentally.   Comparison  of  the  theo- 
retical and  experimental  results  indicates  that  the  theory 
does  predict  the  wave  phenomena  about  the  front  to  a  good 
degree  of  accuracy. 

Conclusions  are  drawn  regarding  the  transient  wave 
response  to  a  wavemaker  and  compared  to  those  of  John  W. 
Miles  regarding  the  transient  gravity  wave  response  to  an 
oscillating  pressure.   The  transient  features  of  the  wave 
front  considered  are  the  width  of  the  transient  zone  at  a 
given  time  since  the  disturbance  to  the  free  surface  of  the 
fluid  and  the  length  of  the  transient  period  at  a  given 
distance  from  the  disturbance.   It  is  shown  that  the  con- 
clusions regarding  transient  gravity  wave  phenomena  apply 
without  qualification  to  the  waves  generated  by  both  oscil- 
lating pressure  sources  and  moving  surfaces. 
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NOMENCLATURE 


6-  Dirac  delta  function 

n-  Wave  height 

g  -  Acceleration  due  to  gravity 

G  -  Green's  function 

r-  Gamma  function 

h  -  depth  of  fluid 

H  -  Heaviside  step  function 
I   -  Imaginary  part  of  a  complex  variable 
R   -  Real  part  of  a  complex  variable 

s  -  Displacement  of  wavemaker  paddle 

S  -  Amplitude  of  displacement  of  wavemaker 
paddle  at  free  surface 

V  -  The  ratio  of  distance  to  time  (x/t) 

v  -  Same  as  V  above. 

w  -  Frequency  of  oscillation  of  wavemaker  paddle 


INTRODUCTION 

The  object  of  this  study  is  to  develop  an  expression 
wnicn  describes  the  local  wave  effects  near  a  wave  front 
advancing  into  calm  water  that  is  uniformly  valid  about  tne 
front.   Limiting  forms  of  this  expression  as  t  -*■  °°  or  x  ->  °° 

were  developed  historically  by  Lamb   and  by  Cauchy  and 

2 

Poisson   respectively  where  t  represents  elapsed  time  after 

an  initial  disturbance  to  the  free  surface  of  the  fluid  and 
x  the  distance  from  the  disturbance.   Lamb's  result  de- 
scribes the  waves  resulting  from  a  steadily  oscillating 
disturbance  where  the  elapsed  time  since  the  disturbance 
was  initiated  is  assumed  to  be  so  large  that  all  transient 
features  have  disappeared  and  only  the  so-called  steady 
state  phenomena  remains.   The  Cauchy-Poisson  expression 
describes  the  initial  propagation  of  waves  out  from  a  dis- 
turbance but  does  not  account  for  the  time  dependent  aspects 
of  tne  disturbance.   The  exact  formulation  of  these  expres- 
sions depends  on  the  nature  of  the  disturbance,  e.g., 
pressure  variations  or  moving  surface  -  and  are  widely 
available  in  the  literature. 


Lamb,  H.,  "On  Deep  itfater  Waves'",  Proc.     Lond.    Math.    Soc.3 
(2) ,  2,  371,  (1904) . 

Cauchy,  A.,  Poisson,  S.  D.   "Memoire  sur  al  theorie  des 
ondes",  Mem.    do    I' Acad.    Roy.    des    Sciences  ,  i.  (1816). 


Neither  of  the  above  theories  attempt  to  describe  the 
transient  wave  phenomena  spanning  their  respective  limits. 
This  study  develops  a  tneory  that  predicts  the  waves  in 
this  transient  zone  and  which  upon  passage  to  the  appro- 
priate limit  yields  the  nistorical  result.   The  waves  gen- 
erated by  an  oscillating  vertical  boundary  (a  wavemaker) 
and  propagating  out  into  an  infinitely  long  channel  are 
investigated  both  theoretically  and  experimentally.   The 
channel  is  assumed  to  be  infinitely  wide  so  that  the 
problem  is  two-dimensional. 

During  the  problem  formulation  phase  of  this  analysis, 

tne  attention  of  the  author  was  directed  to  a  paper  by 

3  . 

Miles  f  which  considered  the  transient  gravity  wave  response 

to  an  oscillating  pressure.   Although  the  mechanism  by  which 

transient  waves  are  generated  is  different  in  Miles' 

approach,  the  qualitative  features  of  the  resulting  waves 

siioula  be  similar  to  tnose  generated  by  a  wavemaker.   In  nis 

paper  Miles  developed  an  expression  predicting  the  transient 

waves  thac  was  uniformly  valid  about  the  wavefront.   His 

solution  method,  though  apparently  consistent,  was  not 

immediately  satisfying  as  it  involved  finding  the  asymptotic 

expansion  of  a  double  integral  by  applying  the  method  of 

stationary  phase  twice  in  succession.   It  was  felt  that  a 


JMiles,  J.  W.,  "Transient  Gravity  Wave  Response  to  an  Oscil- 
lating Pressure",  J.     Fluid   Meek.,    13,  (1962),  pp.  145-150. 
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more  direct  means  of  obtaining  a  uniformly  valid  expres- 
sion was  desirable  -  both  to  verify  Miles '  result  and  to 
obtain  clearer  insight  into  tne  transient  problem. 

This  study  was  undertaken  primarily  in  the  interest  of 
academic  completeness;  viz.,  to  bridge  the  gap  between  the 
solutions  of  Lamb  and  Cauchy-Poisson.   However,  several 
engineering  applications  can  be  foreseen.   The  utilization 
of  wavemakers  in  towing  tanks  with  corresponding  applica- 
tions such  as  ship  model  and  ocean  structure  testing  is  of 
direct  interest  to  the  naval  architect  or  ocean  engineer. 
Transient  features  of  generated  waves  may  allow  some 
aspects  of  model  testing  to  be  explored  more  rigorously  and 
completely.   Further  application  is  found  in  the  area  of 
ocean  wave  forecasting.   The  wave  generating  models  used 
presently  in  predicting  the  wave  phenomena  associated  with 
a  storm  at  sea  are  not  completely  satisfactory  and  a  new 
model,  based  partially  or  entirely  on  the  expression 
derived  in  this  study  could  be  formulated. 
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THEORETICAL  ANALYSIS 

Consider  the  waves  generated  by  a  vertical  plate  in 
oscillatory  motion  given  by  (1)  in  a  fluid  bounded  by  a 
free  surface,  the  vertical  plate,  and  of  infinite  extent 
otherwise.   This  so-called  wavemaker  problem  is  solved  in 
Appendix  3  and  the  resulting  expression  for  the  waveheight 
as  a  function  of  distance  and  time  is  given  by  (2) . 


(1) 


=  S  eW  y/g  sin(wt)  H(t) 


where   S  =  amplitude  of  displacement  at  the 
free  surface 


(2)    n(x,t)  = 


H(t)  is  the  Heaviside  step  function 
wS 


7T 


cos  (ux) 
u  +  w2/g 


sin(wt)  -  sin(/gu  t) 
w  -  /gu 


sin(wt)  +  sin(/gu  t) 
w  +  /gu 


du 


Introducing  the  change  in  variable  u  =  a2/g 


(3)    n(x,t)    =   2WS    p   ocos[(o2/g)x] 


a2    +  w2 


sin(wt)    -    sin(at) 


w    -   a 


sin(wt)    +    sin(at) 
w   +   a 


do  = 


( 


=1 


m 


2  2  2  2 

i(-x+at)         i(-x-at)  i(-x+wt)  i(°-x--wt) 


2wS     a2e 


-^e      g 


-wae 


+wae 


da 


(a   +wz)     (a+w)     (a-w) 


Using  nondimensional  variables  x  =  w2x/g  ,   t  =  wt, 
a  =  a/w  and  n  =  -~   ,  (3)  reduces  to 
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(4) 


n(x,t)   = 


m 


f~2ei (a  2x+at) _a  2ei (a  2x-at) _aei (a2x+t) +aei (a  2x-t) 

(a2+l)  (o+l)     (o-l) 


da 


Equation  (4)  for  n(x,t)  may  be  considered  as  the  sum 
of  four  integrals. 


(5) 


n(x,t)  =  im  {ix  -  i2  -  i  +  i4} 


where 


I,  = 


L  = 


I->  = 


I,  = 


f(a)  e 


i (a2x+at) 


da 


£ ,    x   i (a  x-at)  , 
f(a)  e         da 


r|jo)_ei(a2x+t)do 


rUSL«i  <»**-«  do 


and   f(a)  = 


a*-l 


The  asymptotic  behavior  of  these  integrals  for 
large  x  or  t  is  readily  evaluated  by  treating  a  as  a  com- 
plex variable  and  employing  tne  residue  theorem  and 
Cauchy ' s  integral  theorem  of  complex  variable  theory. 
Their  behavior  for  small  x  or  t  is  given  in  Appendix  H. 
In  the  complex  a    plane,  the  real  integrals  in  (5) 
are  transformed  to  complex  integrals  which  are  to  be 
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evaluated  along  the  contour  consisting  of  the  ray  from  the 
origin  to  infinity  along  the  positive  real  axis. 

It  is  easily  shown  (see  Appendix  C)  that  the  asymptotic 
behavior  of  integrals  I, ,  I~  and  I.    is  given  by  the  follow- 
ing expressions  for  large  t  or  x: 

1±   %   i*  ei(x+t)  +  0(t"3) 

I3  „  ^Lei(x+t)  +  0(t"2) 
i4  *  ^Lei(x-t)  +  Q(t.2) 

The  asymptotic  behavior  of  (5)  then  reduces  to  the 
following: 

(6)  ,MM»j^M-I2 

The  asymptotic  behavior  of  I?  is  not  as  easily  as- 
certained as  was  that  of  I,,  I~  and  I..   In  the  evaluation 
of  these  integrals,  the  best  choice  of  an  equivalent  con- 
tour to  facilitate  the  determination  of  their  asymptotic 
behavior  was  obvious.   If  a  similar  contour  is  chosen  for 
tne  evaluation  of  I2,    however,  the  integral  diverges  along 
portions  of  the  contour.   Consider  then  the  exponential 
term  of  I~  in  expanded  form  where  the  substitution 

i  A 

a  =  r  e    has  been  made.   It  takes  the  form 
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n)  i(a2x-at)     it{r2 (x/t)cos(26) -r  sin(6)}  x 


x  e-t{r2 (x/t)sin(26)-r  sin(9)} 


A  good  choice  for  an  equivalent  contour  would  be  one  on 
which  this  term  becomes  small  as  x  and  t  become  large  and 
vanishes  when  the  contour  is  extended  to  infinity.   For 
this  term  to  become  exponentially  small  for  large  values  of 
t  (or  x) ,  the  following  restriction  on  r  and  0  applies: 

(8)  r2(x/t)  sin(29)  -  r  sin(6)  >  0 

Let  x/t  =  V:   then  for  9  negative   r  <  1/2V  cos (9) 

for  9  positive   r  >  1/2V  cos (9) 
for  9=0       r  =  1/2V 

Thus,  any  contour  with  the  general  features  of  C~ 
shown  below  is  a  satisfactory  equivalent  contour  for  the 
path  from  0  to  +  °°  along  the  real  axis .   The  point 
a  =  1/2V  +  Oi  where  the  contour  C~  crosses  the  real  axis  is 
a  stationary  point  of  the  integrand  and  must  correspond 
to  a  saddle  point  as  the  integrand  is  a  function  of  the 
complex  variable  a  and  therefore  satisfies  Laplace's 
equation. 
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a-plane 


ia 


Here  singularities  of  tne  integrand  of  I~  have  been 
indicated  with  x  marks  and  in  all  cases  are  simple  poles. 
The  contour  along  the  real  axis  is  denoted  by  C,  and  has 
been  indented  at  a  =  1 .   Thus /  the  above  integral  is 
defined  in  the  Cauchy  principle  value  sense.   If  the 
integral  exists  in  the  conventional  sense/  then  its  value 
will  necessarily  agree  with  the  Cauchy  principle  value  so 
there  is  no  loss  in  generality. 

The  existence  of  a  saddle  point  at  a  =  1/2V  suggests 
the  most  convenient  equivalent  contour  on  which  to  evaluate 
I~.   If  C~  is  deformed  so  that  it  coincides  with  the 
steepest  descents  path  from  the  saddle  point/  then  divergent 


integrals  are  avoided  as  the  integrand  of  I-  becomes 
exponentially  small  on  such  a  path  (granting  suitable 
restrictions  on  f(a)).   This  contour  is  shown  below. 


16 


a-plane 
V  <  1/2 


10 


-1 


a-plane 
V  >  1/2 


17 
Portions  of  C~  that  join  the  steepest  descents  path  to 
the  origin  and  to  the  positive  real  axis  at  infinity  have 
been  labeled  C-,  and  CR.   The  location  of  the  saddle  point 
relative  to  the  simple  pole  at  a  =  1/2V  depends  on  the 
parameter  V  =  x/t.   Therefore,  two  general  cases  are  shown. 

The  closing  contours  C,  and  C_.  are  chosen  in  the  fol- 

i.         xv 

lowing  way : 

1)  C-.  is  the  steepest  descents  path  from  the  origin  to 
infinity  and  intersects  the  contour  C~  at  infinity. 

i  ft 

2)  C_.  is  the  arc  of  the  circle  Re    joining  the  contour 

C~    to  the  positive  real  axis. 

Consider  the  case  V  =  x/t  <  1/2.   Application  of 
Cauchy ' s  integral  theorem  and  the  residue  theorem  lead  to 
the  following  expression: 


(9)    I9    -    iTTRes(l)    =    /        +   /        +    f         -    27TiRes(l)    -    2TTiRes(-i) 
Z  Ll  U2  UR 

R+00 


The  asymptotic  properties  of  each  term  in  the  above 
equality  can  now  be  evaluated. 


„   /  • \    lim 
Res(-i)  = 

a-*-- 1 


(g+i)o2   it(o2V-a) 


a4  -1 


i   -t  -itV 
=  j   e   e 


This  term  is  exponentially  small  and  can  be  neglected 
in  comparison  to  terms  algebraically  small. 
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Res(l)  = 


lim 
a-H-1 


(a-l)o2   it(a2V-a) 


a*-l 


i(x-t) 


/r  ->  0  as  R  + 
UR 


It  is  easily  shown  (see  Appendix  C)  that  /„   =  0 (t   ) 

Ul 
and  is  therefore  small  compared  with  other  contributions  to 

I«.   Therefore,  (9)  reduces  to 


(10) 


12   = 


^ei(x-t)  +  0(t"3) 


Consider 


(11) 


it(a2V-a) 


C2  a*-l 


da 


where  V  =  x/t  is  a  parameter  to  be  specified  and  is  less 

than  1/2  for  the  case  under  consideration.   Introducing  the 

1  ? 

change  in  variable  -  ttt  +  it  =  a  V  -a,  the  following  con- 
siderations apply: 


"  w  +  iT  =  j2y  '  a 


a  =  ^  [  1  +  /4lV  t1/2  ] 


~  [  1  +  yi  /2  ]  ,  y  =  /HV 


do 
dx 


+  Y 


4Vt 


71 
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On  C~  in  the  upper  half  plane  and  in  the  half 
plane  to  the  right  of  C~ 


a  =  £    [  1  +  YT  /2  ]  =  a+  , 


Y 


di    4Vt 


75 


On  C2  in  the  lower  half  plane  and  in  the  half  plane 
to  the  left  of  C„ 


=  4r  [  1  "  YT1/a  ]  = 


da 


-Y 


2V 


dT    4Vx 


/2 


Therefore,  the  integral  along  C~    is  given  by 


(12) 


it{a2 (x/t)-a> 


da  = 


C2  a*-l 


it(a2V-a) 


da 


C2  a*-l 


=  e 


-it 
4V 


da  , 

+   -ti  , 
e    di  + 


o  a,-l  di 


da 


-tT, 

e    d  t 


a1*-!  dx 


If  the  Laurent  expansion  about  t  =  0  is  obtained,  the 
above  integrals  may  be  evaluated  by  integrating  term  by 
term.   It  is  more  convenient,  however,  to  subtract  off  the 
pole  at  o  =  1  and  to  then  expand  the  remainder  of  the  ex- 
pression about  i  =  0  so  that  the  effect  of  the  pole  may  be 
considered  separately  with  a  view  towards  obtaining  a  uni- 
formly valid  result.   These  expansions  are  developed  in 
Appendix  D  and  are  given  below. 
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(13) 


da+         -y/8 


a*-l     dx  1-2V 


T*   + 


1-2V 


Y 


+ 


-Vz 


1-2V 


(-)  U-4V-4V2}t     /2         y{  (1+4V2)  2-8V(l+2V)  2} 
(1+2V) (1+4V2)  (1+2V) 2(1+4V2) 2 


(-)  y2{  (1-I-4V2)  3-4V(l+2V)  3  (3-4V2)  }t  /2, 

_  n      n   i  i  '  1  •     •     •      • 


(1+2V) 3     (1+4V2) 


2^    3 


(14) 


a2      da_        y/| 


o*-l     di         1-2V 


Vz 


1-2V 


-1 


-T 


-y2 


1-2V 


+       {1-4V-4V2}  -Vz     .     y{  (1+4V2)  2-8V(l+2V)  2>    + 

(1+2V) (1+4V2)  (1+2V) 2 (1+4V2) 2 


+    Y2(  (1+4V2)  3-4V(l+2V)  3  (3-4V2)  }       l/2 
(1+2V) 3     (1+4V2) 


Noting  that  terms  containing  integral  powers  of  x 
cancel  identically  upon  substitution  into  (12) ,  the 
resulting  expression  becomes 
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(15) 


it(a2V-a) 


da   = 


-y/8    -it/4V 


C2   a*-l 


1-2V 


e        dx 
-r 


o    x 


e        dx 
V2         1-2V 
Y    " 


+    e 


-it    r 
4V 


2|      ^/8    T 
o1^      1-2V 


^2 


/2.    1-2V 

O     T   '    +     

Y 


Y/3[1-4V-4V2] 
(1+2V) (1+4V2) 


_    Y3/8  [  d+4V2)  3-4V(l+2V)  3  (3-4V2)  ]     ^/2    +    0(t3/2)1    e"tTd. 
(1+2V) 3 (1+4V2) 3  J 


The  first  bracketed  term  in  (15)  includes  the  effect 
of  the  pole  at  a  =  1.   The  second  bracketed  term  can  be 
expressed  as  a  series  of  Gamma  functions  and  is  not  in- 
fluenced by  the  pole  at  a  =  1.   Introducing  the  change  of 

V- 

variable  t    =  s,  the  integrals  in  the  first  bracketed  term 
above  are  transformed  to  the  following: 


(16) 


e    dx 
-n i — ^tt 


0   T 


y2+ 


1-2V 

Y 


-tT, 

e    dx 


0   T 


Y, 


1-2V 
Y 


°°      -tsi    ,           r°°      ~ts    -    -j 
e        2sds  e 2sds 

o    .  1-2V  '   I        1-2V 


s  + 


Y 


o  s  - 


Y 


Let  s  =  -s  in  the  second  integral  of  (16) . 

Then 
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(17) 


f~  -ts2 


2sds 


os  +  1-2V 


-ts 


2sds 


s  - 


Y 


1-2V 
Y 


-ts 


2sds 


'°  -ts2 


s  + 


1-2V 
Y 


2sds 


J-c 


s  + 


1-2V 
Y 


-ts 


2sds 


J-c 


s  + 


1-2V 
Y 


=  2 


r°°    ,      1-2V    l-2Vv    .   2 

(s  +  -  )  -tsz-, 

y  y   e ds^ 


J-c 


s  + 


1-2V 

Y 


=  2 


-ts 


ds  - 


1-2V 
Y 


-ts 


ds 


I    s  + 

J—  oo 


1-2V 
Y   -1 


y2 

variable  y  =  t 


With  the  change  in  variable  y  =  t/zs  the  first  inte- 
gral in  (17)  is  readily  evaluated  and  the  second  placed 
in  a  more  convenient  form.   Equation  (17)  then  reduces  to 


(18) 


2  f/¥  fV2  +  ?  (±=2L)    W(z)  1 

L  l     Y  J 


r°° 


where   W(z)  =  — 


rv 


z  -  y 

'— oo  J 


[I  (z)  >  0] 

m 


=  _tV2  (i±) 

Y 


W(z)  can  be  represented  by  the  following  expression  : 


-z 


W(z)  =  e    erfc(-iz) 


where   erf c  (z)  = 


/? 


e  fc  dt 


Abramowitz,  M.  and  Stegun,  I.  A.  ed.,  Handbook    of  Math. 
Functions ,    New  York,  (1965) ,  p.  297. 
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On  substituting  for  W(z) ,  (18)  can  be  placed  in  the 


form 


(20) 


/i  t"1/2  +  I  (±=2V,  w(z) 

i  y 


=    2 


/it    t    '  *■    -    i  it  ( )  e 


-<±^Vt 


er 


f0  {i(iz2v,tlA} 


Therefore,  (11)  can  be  expressed  as  follows  for  the 
case  V  =  x/t  <  1/2: 


(21) 


'    a2    it(a2V-a) ,   A 
e         da  'v 

C2  a1*-! 


-it 


ZlZi  e  4V 
1-2V 


2/^t~/2-  2TTi(i^V)e    Y 
Y 


-<i^)2t 


er 


fc{i(irp)t^} 


+  e 


-it/4V 


y>* 


Y/4    ,1,   -l/2         y/4r(l/2)  {l-4V-4V2}t  /2  , 
1-2V  1  K2}     C  (1+2V) (1+4V) 


+  -Y3/4{  (1+4V2)  3-4V(l+2V)  3(3-4V2)  }T(3/2)t  /z  +  0(t-5/2) 

(1+2V) 3 (1+4V2) 3 


Substitution  of  (21)  into  (10)  leads  to  the  following 
expression  for  (6) : 
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(22) 


,  ±7T        i  (x-t) 

n(x,t)    *  Im  .-j-  e 


-it       ,1-2V, 2^ 

i?T         4V       ^     V 
-j-  e        e 


erfc{i(i^)t1/2}  + 


+    e 


-i(k-?»  rv{i-4v-4v2} 


(1+2V) (1+4V2)  r( 


1/2)  t  /a  + 


+  2ivV2{(l+4VV-4V(U2V)3(3-4V^)}r(3)t-V2+  0(^)1}  <V<i) 
(1+2V) 3(1+4V2) 3         ^  J)  z 

where  y   =  /4iV  =  V4i  x/t 

Consider  the  case  V  =  x/t  >  1/2.   Application  of 
Cauchy's  integral  theorem  and  the  residue  theorem  leads  to 


(23) 


I«  -  iff  Res(l)  = 


+ 
°1 


2  "R 

R->oo 


Similar  considerations  regarding     and 

Ci 


'R 


apply  as  for  the  case  V  =  x/t  <  1/2.   Therefore,  the 
asymptotic  behavior  of  L  for  large  t  or  x  is  given  by 


(24) 


1^  = 


it(a2V-a) 


c2  a*-l 


da  +  -r- 


X7T   i  (x-t)  ,  rt  ,  -3 


+0(t  J) 


Proceeding  as  for  the  case  V  <  1/2,  introduce  the 
change  in  variable  -  j^   +  ix  =  a2V  -  a .   All  previous  con- 
clusions regarding  the  branches  and  branch  cuts  of o(x) 
apply,  hence 
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(25) 


it (a2V-a)  _ 
e         da  = 


C2  0»-l 


=  e 


■it 
4V 


2 
+ 


■"o  o+-l     dx 


+   -tx  , 
e    dx  + 


oo    g 


da 


-  -tx, 

—  e   dx 


o  a  -1  dx 


The  Laurent  expansions  previously  obtained  about  x  =  0 
are  valid,  hence 


(26) 


\«  e 


-it 
4V 


Y/4  r(l/2) 


t"^ 


Y/4{l-4V-4V2}r(l/2) 


t-^ 


1-2V 


(1+2V) (1+4V2) 


Y  3/4{  (1+4V2)  3-4V(l+2V)  3  (3-4V2)  }  r{3/2)t~ '%  +  0(t  /z) 
(1+2V) 3 (1+4V2) 3 


-it 

y/8     4  V 
_iz e 

1-2V 


oo    -ts 


2sds 


s    1-2V 

o     +  


oo    -tS 


2sds 


1-2V 

Y 


The  remaining  integrals  are  treated  in  a  similar  manner 
as  before  except  that  the  substitution  s  =  -s  should  be  made 
in  the  first  integral  in  lieu  of  the  second  integral  of  (26) 
so  that  the  resulting  formulation  exists. 


(27) 


-ts 


sds 


00      -ts*     , 
e  sds 


0    s    + 


1-2V         j 
Y 


0    s    - 


1-2V 

Y       -1 


=    2 


-ts2     , 
e  sds 


_oo     s     - 


1-2V 
Y 


sds 


0     s    - 


1-2V 
Y       J 


=    2 


e  sds 


'-co    s      - 


1-2V 
Y 


=    2 


-ts 


ds   + 


J  -c 


1-2V 
Y 


-ts 


ds 


J- 


00     s 


1-2V 
Y       -1 


1/ 
Let   y   -    t  /2s. 


=    2 


?  t-V2_^    (lz2V 


/? 


(±-J^L)      W(Z) 

1  Y 


i/2/l-2V) 

where      z   =   t/z( 

Y 


(28) 


=    2 


/1-2V,  2. 
/tt   t     /z   +    itt( )e  '  erfc{-i(— — )t'    }     1 
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The  evaluation  of  (25)  is  now  complete  and  the  wave- 
height  n(x,t)  is  therefore  given  by 
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(29) 


n(x,t)   ^  i 


m 


-it2       (1-2V    2 

— r-  e  e  '  erfc{-i( )t/    }    + 

4  Y 


+    e 


4~x         4 


rJ: 


V/2/2{l-4V-4V2}r(l/2)     t-y2    + 
(1+2V) (1+4V2) 


,     2iV  ^2{  (1+4V2)  3-4V(l+2V)  3(3-4V2)}    p.3>     -3/2 
(1+2V) 3     (1+4V2)3  ^ 


+  o(t  5/2 


(t     ") 


(V    >    1/2) 


where     y    =    ^4iV 


The  asymptotic  expansions  for  n(x,t)  have  been  derived 
for  the  two  separate  cases  of  V  =  x/t  <  1/2  and  V  =  x/t >l/2 
That  these  two  expansions  are  equivalent  can  be  shown  by 
exploiting  either  of  the  following  symmetry  relations: 


(30) 


erfc(-z)  =  2  -  erfc(z) 


-z 


(31)  w(-z)  =  2e     -  w(z) 

where   w(z)  =  i/rr 


e    /(z-t)dt   (I   z  >  0) 

J   —  no 


Considering  equations  (22)  and  (29)  for  ri(x,t), 
symmetry  relation  (30)  is  most  convenient. 

From    (30) , 
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(32) 


er 


,     r  .  ,1-2V»  ,  Vzx         n  *    r     .  /l-2Vva.y2i 

fc{i( )t/   }    =    2    -   erfc{-i( -)  t /    } 


Substituting    (32)    into    (22)    one    obtains 


(33) 


iTT      i  (x-t) 

n(x't)   Mi      Te 


ITT  4x  V 

-r-  e  e  ' 


-i<i=2)  t'/2 


2    -   erfc{-i( )    t      i 


-i( 


+   e 


4x 


-) 
4; 


t-v2 


0(t     ") 


(V   <    1/2) 


Therefore 


/      *.i     ,    t      /    i^         4x         l     y  *    r     •  #l-2Vv       */2*.l. 

n(x,t)    %    Iffl        j  e  e  T  erfc{-i(  — — )    t/2>  + 


+    e 


•  tt  tTn 

"1(4^   "    4)     r 


0(t    l/l)  \         (V   <    1/2) 


However,  (34)  is  identical  to  (29)  -  the  case  V  >  1/2. 
Therefore,  the  restriction  V  <  1/2  is  unnecessary  and  can 
be  removed  and  either  expression  completely  describes  the 
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waves.   The  final  result  then  is  uniformly  valid  about 
the  wavefront  whose  position  is  defined  by  V  =  x/t  =  1/2. 
Some  simplification  of  n (x  t)  is  possible  and  there  are 
many  equivalent  expressions  for  it.   A  particularly  compact 
expression  is  given  below. 


ITT 
(35)  ,  .    v      „      T  flTT         i(X-t)  jrr^-/T  /->      \  4       n 

n(x,t)    ^   I      <    -—  e  erfc [/x (l-t/2x) e        ]    + 

"i(4^        4*      ~  /ix7t{l-4x/t-4x2/t2)t"1/2     t     Q(t-3/2)" 

2(l+2x/t)     (l+4x2/t2) 


For  the  cases  V  <<  1/2,  V  >>  1/2,  and  V  =  1/2,  expres- 
sions for  the  waveheight  n(x,t)  can  be  developed  in  an  al- 
ternative manner.   The  effect  of  the  simple  pole  at  a  =  1 
is  either  negligible  or  an  easily  evaluated  residue  contri- 
bution in  these  cases.   These  alternative  expansions  for 
n(x,t)  have  intrinsic  merit  and  in  addition  provide  a  check 
on  the  validity  of  the  uniformly  valid  expression  that  has 
been  derived.   They  are  available  in  Appendix  E. 
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EXPERIMENTAL  PROCEDURE 

Experiments  were  conducted  in  the  Ship  Model  Towing 
Tank  located  in  Building  48,  Hydrodynamics  Laboratory.   The 
purpose  of  the  experiments  was  to  provide  a  comparison  be- 
tween the  waves  generated  by  a  wavemaker  and  those  pre- 
dicted by  the  derived  theoretical  expression. 

Control  of  the  oscillatory  motions  of  the  wavemaker 
paddle  was  accomplished  by  utilizing  the  sinusoid  voltage 
output  of  a  function  generator  as  the  input  to  the  paddle 
motion  control  mechanism  -  a  hydraulic  device.   A  voltage 
output  proportional  to  the  deflection  of  the  wavemaker 
paddle  was  fed  back  and  compared  to  the  input  voltage  pro- 
viding close  synchronization  of  the  paddle  motion  to  the 
desired  motion.   The  amplitude  of  horizontal  displacement 
oi  the  paddle  at  the  water  surface  was  measured.   Both 
frequency  of  paddle  oscillations  and  maximum  amplitude  of 
displacement  were  controllable  parameters  throughout  the 
experiment.   It  was  not  found  necessary  to  vary  the  maximum 
amplitude  of  displacement  from  its  initial  value  of  1  inch, 
however.   Paddle  frequency  was  varied  from  0.8  cycles  per 
second  (cps)  to  2.0  cps .   Above  a  frequency  of  1.2  cps  non- 
linear effects  were  encountered  and  severe  vibrations  of 
the  wavemaker  machinery  noted.   Transverse  wave  patterns 
were  also  induced  indicating  that  significant  three- 
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dimensional  features  were  present  in  the  flow.   As  a 
result,  all  experimental  data  for  runs  above  1.2  cps  was 
not  evaluated.   Experiments  were  not  conducted  for  fre- 
quencies below  0.8  cps  so  as  to  avoid  the  effects  of 
finite  towing  tank  depth. 

The  height  of  the  waves  generated  by  the  wavemaker 
was  measured  with  waveprobes .   These  instruments  are  var- 
iable resistance  devices  whose  electrical  resistance  changes 
with  depth  of  immersion  into  the  water.   The  voltage  output 
of  the  waveprobes  was  connected  to  the  input  of  a  strip 
chart  recorder  (Sanborn,  series  7700) .   Thus  a  continuous 
record  of  the  waves  at  a  particular  location  from  the  wave- 
maker  paddle  versus  time  is  obtained.   Calibration  of  the 
waveprobes  was  accomplished  by  immersing  the  probes  a 
measured  distance  into  the  water  and  adjusting  the  deflec- 
.tion  of  the  strip  chart  recorder  scribe  to  a  convenient 
value.   The  deflection  of  the  wavemaker  boundary  was  also 
converted  to  a  voltage  signal  and  connected  to  the  chart 
recorder  input  so  that  the  motions  of  the  wavemaker  were 
continuously  monitored. 
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RESULTS 

To  facilitate  comparison  of  the  waves  predicted  by 
equation  (35)  and  those  measured  experimentally,  all  dimen- 
sional quantities  are  normalized  to  a  consistent  base  and 
the  results  illustrated  graphically.   The  transformation 
equations  relating  dimensional  to  nondimensional  quantities 
are  those  used  in  the  theoretical  analysis  in  the  derivation 
of  equation  (35)  and  are  as  follows: 

2  — 
W  X 

x  =  where  w  =  frequency  of  oscillation 

g  of  wavemaker 

t  =  w  t  g  =  acceleration  due  to  gravity 

—  S  =  amplitude  of  oscillation  of 

n  =  ~— ^r  wavemaker  at  the  free  surface 


2  S 


All  bar  quantities  are  dimensional 


Predicted  values  of  the  wave  height  are  readily  com- 
puted manually  using  tables  to  determine  values  of  the  error 
function  complement  for  complex  arguments.   Equation  (35)  is 
well  suited,  however,  for  computer  evaluation  of  r\(x,t). 
Computer  programs  for  this  purpose  are  presented  in  Appendix 
G.   In  addition,  a  function  subprogram  for  the  error  func- 
tion complement  for  complex  arguments  is  listed  separately. 

Computer  facilities  were  also  used  to  plot  the  theo- 
retical and  experimental  values  of  wave  height  using  the 


11 

Stromberg-Carlson  4020  (SC-4020)   cathode  ray  tube  plotting 
device.   Graphs  of  experimental  results  are  computer  faired 
with  points  corresponding  to  normalized,  measured  values 
of  wave  height  encircled.   Comparison  of  faired  to  experi- 
mental results  show  that  the  faired  results  differ  very 
little  from  measured  values.   Reproductions  of  the  wave 
records  obtained  during  the  experiments  are  available  in 
Appendix  F. 

Figures  1  through  5  compare  theoretical  and  experimental 
results  for  five  values  of  x.   Wave  height  n(x/t)  is  plotted 
versus  time  t  for  fixed  values  of  distance  x.   These  figures 
therefore  show  the  timewise  distribution  of  wave  height  as 
seen  by  an  observer  located  at  a  fixed  distance  from  the 
wavemaker.   Figure  6  is  generated  from  equation  (35)  for 
x  =  65  and  illustrates  well  all  the  prominant  features  of  the 
timewise  distribution  of  waves.   Figure  7  shows  the  spacewise 
distribution  of  wave  height  for  various  of  time  t.   These 
figures  are  arranged  sequentially  and  show  the  wave  front 
propagating  out  from  the  wavemaker  with  increasing  time. 
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Fig.  1-A.   Wave  Height  vs.  Time  for  x  =  8.7687.   Based 
on  Equation  ( 35) . 
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Fig.  1-B.   Wave  Height  vs . » Time_f or  x  =  8.7  68  7 

on  Measured  Values  (x  =  11.6667  ft, 
f  =  0.  7826  cps)  . 
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ig.  2-A.   Wave  Height  vs.  Time  for  x  -  13.6493.   Based 
on  Equation  (35). 
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Tig.  2-B.   Wave  Height  vs.  Time  for  x  =  13.6493 
on  measured  values  (x  =  11.6667  ft, 
f  =  0  .9754  cps)  . 
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Fig.   3-A.   Wave  Height  vs.  Time  for  x  =  19.7632.   Based 
on  Equation  (35). 
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Fig.  3-r>.   Wave  Height  vs.  Time_for  x  =  19.7632.   Based 
on  measured  values  (x  =  11.6667  ft, 
f  =  1. 1749  cps)  . 
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Fig.  4-A.   Wave  Height  vs.  Time  for  x  =  29.4239.   Based 
on  Equation  (35). 
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Wave  Height  vs.  Tim e_ for  x  =  29.4239 
on  measured  values  (x  =  38.1667  ft, 
f  =  0. 7926  cps) . 
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Fig.  5-A.   Wave  Height  vs.  Time  for  x  =  44.8634.   Based 
on  Equation  (35 ) , 
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Fig.  5-B.   Wave  Height  vs.  Time_for  x  =  44.8634.   Based 
on  measured  values  (x  =  38.1667  ft, 
f  =  0.9787  cps) . 
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Fig.  6 


Wave  Height  vs.  Time  for  x  =  65.   Based  on 
Equation  ( 35) . 


Figure  6  illustrates  well  all  the  prominant  features  of 
the  timewise  distribution  of  waves.   The  wave  height  in- 
creases monotonically  to  a  maximum  value  and  then  enters  an 
oscillatory  phase  until  the  steady  state  wave  amplitude  is 
reached.  ! 

Figure  7  shows  the  spacewise  distribution  of  wave 
height  for  various  values  of  time  t.   These  figures  are 
arranged  sequentially  and  show  the  wave  front  propagating 
out  into  calm  water  with  increasing  time. 
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Fig.  7-A 


Wave  Height  vs.  Distance  for  t 
Based  on  Equation  (35). 
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Fig.  7-B.   Wave  Height  vs.  Distance  for  t 
Based  on  Equation  (35). 
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Fig.  7-C.   Wave  Height  vs.  Distance  for  t  =  90.0 
Based  on  Equation  (35). 
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Fig.  7-D.   Wave  Heigh-t  vs.  Distance  for  t 
Based  on  Equation  (35)  . 
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Fig.  7-E.   Wave  Height  vs.  Distance  for  t 
Based  on  Equation  (35). 
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Fig.  7-F.   Wave  Height  vs.  Distance  for  t  =  240.0. 
Based  on  Equation  (35). 


DISCUSSION  — 


The  waves  predicted  by  the  equation  (35)  for  the  wave 
height  as  a  function  of  distance  and  time  compare  well  with 
measured  values.   To  an  observer  located  a  fixed  distance 
from  the  wavemaker,  the  wave  height  increases  monotonically 
with  time  to  a  maximum  value  and  then  enters  an  oscillatory 
phase  during  which  the  steady  state  wave  amplitude  is 
approached.   Both  measured  and  predicted  wave  records  show 
this  characteristic  behavior. 

Measured  and  predicted  values  of  the  wave  height  cor- 
respond very  closely  in  phase,  except  for  very  small  values 
of  t  when  equation  (35)  is  no  longer  valid.   However,  meas- 
ured values  of  wave  amplitude  are  consistently  larger  than 
predicted  values.   This  difference  can  be  attributed  to  the 
boundary  condition  imposed  on  the  moving  surface  of  the 
wavemaker  in  the  theoretical  development.   For  convenience 
of  solution,  this  boundary  condition  was  assumed  to  be  the 
following : 

(1)         s  =  S  eW  y/gsin(wt)  H(t)    (—  <  y  <  0) 

However,  a  more  realistic  approximation  to  the  flow 
near  the  towing  tank  wavemaker  paddle  is  given  by 


(36)        s  =  { 


S(y  *  h)  sin(wt)  H(t)     (y  >  -h) 

(y  <  -h) 
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This  is  not  the  exact  boundary  condition  but  is  a  reason- 
able approximation  to  it.  It  assumes  that  displacements  of 
the  wavemaker  from  the  vertical  are  small  with  respect  to 
the  depth  of  the  wavemaker  padale  (approximately  equal  to 
the  depth  h  of  the  towing  tank  wave  generator)  and  is  the 
linear  approximation  to  the  actual  displacement. 


Wavemaker 


Free  Surface 


In  addition,  the  theoretical  expression  for  n(Xrt)  is  for 
the  case  of  infinite  depth  water.   The  waves  generated  by 
the  towing  tank  wavemaker  have  been  chosen  so  that  finite 
depth  effects  have  a  negligible  effect  on  the  wave  features 
Therefore,  the  boundary  condition  given  by  (36)  is  a  simpli- 
fied model  of  the  physical  situation  shown  above. 

A  tneoretical  development  based  on  the  more  realistic 
approximation  to  the  boundary  condition  poses  no  additional 
complication  to  the  problem  but  does  significantly  add  to 
the  algebraic  tedium  and  perhaps  obscures  evaluation  of  the 
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results  somewhat.   Of  primary  concern  is  the  determination 
as  to  just  how  the  use  of  boundary  condition  (1)  above  in 
lieu  of  (36)  has  influenced  the  problem.   Physical  consid- 
erations suggest  that  assuming  an  exponential  decrease  in 
fluid  velocity  with  depth  instead  of  the  more  realistic 
linear  decay  should  result  in  wave  amplitude  predictions  of 
smaller  magnitude  than  measured  values.   As  less  energy  is 
imparted  to  the  fluid  near  the  moving  surface  (average 
amplitude  of  oscillation  being  smaller) ,  wave  amplitudes 
must  be  correspondingly  smaller  as  energy  considerations  are 
ultimately  related  to  the  square  of  wave  amplitude.   This  is 
substantiated  by  experimental  results. 

Partial  verification  of  the  physical  argument  presented 
above  can  be  readily  obtained  in  spite  of  the  algebraic 
complexities  of  the  more  realistic  boundary  condition.   The 
steady  state  wave  amplitude  can  be  easily  obtained  by  let- 
ting t  approach  infinity  and  evaluating  the  remaining  terms. 
For  the  case  at  hand  this  amounts  to  determining  the  residue 
contributions  to  the  integral  representation  for  the  wave 
height. 

The  solution  to  the  unsteady  wavemaker  problem  (in- 
finite depth  water)  can  be  placed  in  a  more  general  form 
where  the  functional  dependence  of  the  wavemaker  paddle 
displacement  on  depth  is  left  unspecified  (subject  to  the 
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restrictions  previously  mentioned  on  this  function) .   One 
formulation  of  it  is  the  following: 

(37) 


n(x,t)  = 


ws 

TT 


cos (ux) 


eUyF(y)dy 


sin  (wt) -sin  (/g~u  t) 
w  -  /gu 


sin (wt) +sin (/gu   t) 

w    +    /gu 

For   F(y) 

w2 
-  e  s    , 

du 


eu*F(y)  = ± 

o  u  +  w2/g 


and  equation  (2)  is  rederived. 


y  +  h 


h 


For  F(y)  =  - 


0 


y  >  -h 


y    -h 


eUyF(y)dy  =  -  -A-  ( 
0  hu2 


uh  +  e"Uh  -  1) 


Therefore,  the  wave  height  n (x,t)  is  given  by 


(38)  n(x,t)  = 


2ws     (uh+e    -1) 

TT    ■'o        hU2 


cos (ux) 


"/gu  sin(/gu  t) -w  sin(wt) 
gu  -  w2 


du 
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Introducing  the  change  in  variable  u  =  Q2/g   and 
using  non-dimensional  variables  x  =  w2/g  x,  t  =  wt,  a  =  a/w, 
h  =  w2/g  h,  and  n  =  Trn/2s  one  obtains: 


(39) 


n(x,t)  =  I 


m 


.00         j 

{a2h+e~a  h-l}  /  „  i(a2x^at) 


haMa+1)  (a-1) 


ae 


aei(a2x-at)  _  ei(a2x+t)  +  ei(a2x-t)  ,  do 


Following  methods  similar  to  those  used  in  deriving 
(35)  ,  as  t  ■*■  °° 


(40)  n(x,t)  ^  Im  {  iTr  Res[l]  } 

.  (h  +  e"h  -  1)  7T     ,   .  v 
^ cos (x-t) 


where  h  =  depth  of  towing  tank. 

For  the  measured  depth  of  43.5  inches  (3.62  feet)  and 
wavemaker  paddle  frequency  of  0.9  787  cps ,  h  =  4.25. 

-4 .  25 

(41)  n(x,t)  *  (4<25  +  e —  x  TTCOs(x-t)  s  2.4  cos(x-t) 

4.25 


Equation  (41)  compares  favorably  with  measured  steady 
state  values  of  wave  amplitude  corresponding  to  a  paddle 
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frequency  of  0.9787  cps .   These  experimentally  obtained 
values  of  wave  height  are  given  by: 

(42)  n(x,t)  -  2.35  cos (x-t) 

The  ratio  of  the  steady  state  wave  height  associated 
with  boundary  condition  (36)  to  that  associated  with 
boundary  condition  (1)  is  2.4/1.57  =  1.53.   If  this  ratio 
is  assumed  to  hold  in  a  very  loose  sense  throughout  the 
transient  zone,  then  the  predicted  values  of  wave  height 
are  very  close  to  measured  values. 

For  large  values  of  t  or  x,  equation  (35)  for  n(x,t) 

must  approach  the  historical  solutions  of  Lamb  and  of 

Cauchy  and  Poisson.   This  property  can  be  shown  by  expanding 

the  erfc  for  large  positive  and  negative  values  of  its  ar- 
gument.  Letting  t  approach  °°  one  obtains: 

(43)  n(x,t)  *   tt/2  cos  (x-t) 

These  are  the  outgoing  waves  at  infinity  that  Lamb  derived. 
Letting  x  approach  oo  one  obtains  as  the  first  order  term: 

h      _  i/ 

(44)  n(x,t)    *      5giy    x    /2    cos(gt2/4x   -    tt/4) 

2wtt  /2 

Equation  (44)  exhibits  the  general  features  of  the 
Cauchy/Poisson  solution. 
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CONCLUSIONS 

In  his  analysis  of  the  transient  gravity  wave  response 
to  an  oscillating  pressure,  Miles   described  many  of  the 
qualitative  features  of  the  generated  waves  as  were  observed 
in  the  towing  tank  and  predicted  by  the  derived  expression 
for  n(x,t),  Equation  (35).   The  wave  amplitudes  associated 
with  a  concentrated  oscillating  pressure  also  show  the  mono- 
tonic  increase  to  a  maximum  value  followed  by  an  oscillatory 
phase  to  the  steady  state  amplitude. 

In  addition,  Miles  defines  a  rise-time  T  as  the  time 
for  the  timewise  envelope  to  rise  from  0.1  of  its  steady 
state  value  to  its  first  maximum,  and  a  width  X  of  the 
transition  zone  for  the  spacewise  envelope  as  the  distance 
between  the  point  where  the  wave  height  reaches  0.1  of  its 
steady  state  value  and  its  maximum  amplitude.   The  expres- 
sions for  T  and  X  are  as  follows: 


(45)  T  *  9  (x/g) 

(46)  X  -  3  g  w"3/2t1/2 


It  is  of  interest  to  determine  if  the  theory  developed  for 
the  transient  wave  response  to  a  wavemaker  supports  these 
observations  or  modifies  them  in  any  meaningful  way. 


Miles,  Op.     Cit.  ,    pp.  149-150. 
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Equation  (35)  gives  the  transient  wave  response  to  a 
wavemaker  and  can  be  rewritten  as  follows: 


(47)  n(x,t)  *  Im  fif  ei(x~t)erfc[/x~(l-t/2x)ei7r/4]  +  0(t  ^2) 


-  y? 

Terms  of  order  0(t    )  can  be  neglected  for  sufficiently 
large  t  (or  x) .   It  is  convenient  to  treat  the  remainder  of 
(4  7)  as  an  outgoing  plane  wave  with  amplitude  and  phase 
given  functions  of  x  and  t.   With  this  view  in  mind,  (47) 
can  be  rewritten  as  follows: 

(48)  n(x,t)    ^  A  cos[(x-t)    +    <J>] 
where 

(49)  A  -    tt/4     I    erfc[/^(l-t/2x)e17T'/4]     | 

(50)  <l>   =   Arg    {    erf  c  [  /x~(l-t/2x)  e17T/4]     } 

Expressions  corresponding  to  (45)  and  (46)  can  be 
developed  for  the  transient  wave  response  to  a  wavemaker 
by  investigating  the  behavior  of  A  and  using  Miles ' 
definitions  of  X  and  T. 

The  behavior  of  A  as  a  function  of  U  =/x(l-t/2x)  is 
shown  below. 
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/x(l-t/2.x) 


Fig.  8.   Timewise  Wave  Envelope 


A  is  equal  to  0.1  of  its  steady  state  value  (A  =  0.15708  = 
0.1* it/2)  when  U  =  +2.79  5.   A  reaches  its  maximum  value  when 
U  =  -1.5255.   With  tnese  values,  the  rise-time  T  and  the 
width  of  the  transient  zone  X  are  readily  determined.   In 
terms  of  x  and  t,  the  maximum  and  0.1  steady  state  ampli- 
tudes occur  when: 


(51)     /x  (l-t/2x)  =  -1.5255   for  A  =  A, 


max 


(52)  /x    (l-t/2x)    =   2.795  for   A  =    O.lxrr/2    =  AQ>1 


ss 
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For   x   =   x    : 
o 


tm   =    2xQ+    (2) (1.5255)xq 


Vz 


t  =    2xQ-    (2) (2.795)x    /z 

O.iss  ° 


Therefore 


T   =    t      -    tn    ,     „   =    8.641x    /2 
m  0 . lss  o 


(53)  T    a    9x 

o 


For    t   =    t 
o 


Vz 


^m'^V^m'  =   -1-5255      for  A  =  Afflax 

^0.1ss(1-V2x0.1ss)    "    2-"5  for  A  "  A0.1ss 


These    relations    for   x      reduce    to 

m 


x2    -    {t   +(-1.5255) 2}x      +    t2/4    =    0 
m  o  mo 


x0.iss   "    {to+(2-795)2)x0.1ss   +   fco/4   =   ° 


x      and   x„    ,         are    then   given  by 
m  O.lss  J  J 


{t   +(-1.5255)2}    -    1.5255/2t   +(-1.5255) 2 
o  o 

X         = 

m  2 

{t   +(2.795)2}    +    2.795/2t   +(2.795) 2 
o o 

x0.1ss 
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Therefore, 


X  =  xA  -,    -  x 

O.iss    m 


5.495  +  2.795/2t  +7.82  +  1.5255/2t  +2.325 

o  o 


For  sufficiently  large  t  : 

(54)  X  =  3.06  t  /z  -  3  t  /z 

o        o 


Equations  (53)  and  (54)  can  be  placed  in  the  same  form 
as  (45)  and  (46)  by  introducing  dimensional  quantities,  thus 
obtaining: 


x   i 

(55)  T  a  9(  — ) 

g 


o.Vi 


(56)  X  *  3gw  /2t  1/? 


The  theory  developed  for  the  transient  wave  response 
to  a  wavemaker  thus  supports  the  observations  of  Miles  with 
respect  to  the  asymptotic  behavior  of  the  wave  envelope 
A(x,t).   The  conclusions  regarding  the  rise-time  T  and  the 
width  of  the  transient  zone  X  apply  without  qualification 
to  the  transient  waves  generated  by  both  oscillating  pres- 
sure sources  and  moving  surfaces. 
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RECOMMENDATIONS 



Several  subjects  for  future  analysis  are  recommended 
which  are  related  to  the  problem  considered  by  this  study. 
They  include  the  following: 

1.  Develop  a  uniformly  valid  expression  for  the 
transient  wave  response  to  a  wavemaker  similar  to  equation 
(35)  ,  but  using  the  more  realistic  approximation  to  the 
flow  near  the  wavemaker  paddle  given  by  equation  (3  6)  in 
lieu  of  (1) .   Compare  the  waves  given  by  this  expression  to 
the  experimental  results  obtained  in  this  study. 

2.  Consider  the  problem  of  the  transient  gravity  wave 
response  to  an  oscillating  pressure  treated  by  Miles  but 
using  the  method  developed  in  this  study  -  first  obtaining 
an  integral  representation  for  the  wave  height  using  Green's 
theorem  and  an  appropriate  Green's  function. 

3.  Consider  the  transient  wave  response  to  a  wave- 
maker  in  a  fluid  with  finite  depth. 

4.  Consider  the  transient  wave  response  to  an  oscillat- 
ing pressure  in  a  fluid  with  finite  depth. 


APPENDIX  A  £! 

THE  TIME-DEPENDENT  GREEN'S  FUNCTION 


The  solution  of  the  general  problem  of  the  infin- 
itesimal, irrotational,  unsteady  motion  of  a  fluid  with  a 
free  surface  cannot  usually  be  given  in  explicit  form; 
however,  it  can  be  shown  that  the  problem  can  be  reduced  to 
that  of  finding  an  appropriate  Green's  Function.   The 
essential  tool  used  in  solving  the  problem  is  the  time- 
dependent  Green's  Function,  G (x,y , z; £,n , £; t,x) ,  which  is  of 
intrinsic  interest  itself.   It  represents  the  velocity 
potential  which  is  the  solution  to  the  following  problem: 
At  time  t  =  t  a  disturbance  in  the  nature  of  a  source  is 
originated  at  a  point  (x,y,z)  which  coincides  with  (£,n,C). 
It  is  assumed  that  the  pressure  on  the  free  surface  is 
always  zero  and  that  no  immersed  surfaces  are  present.   G 
satisfies  homogeneous  boundary  conditions  on  all  boundaries 
and  appropriate  conditions  at  infinity.   After  G  is  con- 
structed, the  velocity  potential  for  the  given  problem  can 
be  represented  in  terms  of  initial  conditions  and  boundary 
conditions.   If  no  submerged  bodies  are  present  an  explicit 

solution  is  obtained;  otherwise  the  problem  is  reduced  to 

g 

solving  an  integral  equation.   The  two-dimensional  time- 
dependent  Green's  Function  will  now  be  developed  for 
infinite  depth  fluid  and  fluid  of  finite  depth  h. 


Finkelstein,  A.  B.,  "The  Initial  Value  Problem  for 
Transient  Water  Waves",  Comm.    on    Pure    and   Applied  Matin.  t 
Vol.  X,  pp.  511-522. 


G(x,y;g,n,;t,T)  -  Infinite  Depth  Fluid.  56 

G  must  satisfy  the  following  governing  equations: 

o2~  32G  32G  r  ,  r..,  v  -oo<y<0 

V  G  =  ^  +  „-  =  6(x-U6(y-n) 

— °o<X<+°° 

(here  6  represents  the  Dirac  Delta  Function) 

G..  +  gG   =0     on  y  =  0,  the  free  surface 

y  boundary  condition 

G  ■>  0  as  y  ■>  -°° 

G  =  G  =  o        for  t  =  t,  y  =  0,  the  initial 

conditions 

Introducing  the  Fourier  transform  f  of  G 

r00 
f  =    Ge    dx 

J  _oo 

one  obtains 

-  r°° 

f  e    du 


G-i- 

2tt 


Then  substituting  for  G  its  Fourier  transform  representa- 
tion, the  governing  equations  become 

1)  elu^  (f    -  u2f)  =  <5(y  -  n) 

2)  f .   +  gf   =  0      on  y  =  0 

tt    3  y  2 

3)  f  •>  0  as  y  -*■  -00 

4)  f  =  f   =  0        for  t  =  t,  y  =  0 

The  general  solution  to  the  above  problem  can  be  ex- 
pressed in  the  following  alternative  ways: 

elu^f  =  aeuy  +  be"Uy  -  a'cosh(uy)  +  b'sinh(uy) 
=  a"cosh[u(y+h) ]  +  b"sinh [u (y+h) ] 
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for  y  <  r\  ,    solve  equation  (1)  subject  to  equation  (3)  only. 

iu^j.       |u|y 
e    f  =  a  e1  IJ 

for  y  >  n /  solve  equation  (1)  subject  to  equation  (2)  only. 

iuE,.        ,  ,   v    ctt    ,  ,   x 
e   *f  =  c  cosh(uy)  -  -^j-  smn(uy) 


Require  f  to  be  continuous  at  y  =  n  and  that  the  derivative 
of  f  have  a  discontinuity  of  magnitude  -1  at  y  =  r\,    that  is 


e   sf 


iuE£ 
-  e   ^f 


y=n+ 


=  -1 


y=n- 


From  the  continuity  requirement  above 


e'u^[c  cosh(un)  -  sinh(un_)]     (y  <  n) 


e   sf  =  i 


gu 


4-+- 

e'  '^[c  cosh(uy)  -  sinh(uy)]     (y  >  n) 


From  the  requirement  on  the  derivative  of  f 


ctt  +  guc  =  g 


subject  to    c  =  c.  =  0    for  t  =  t 


The  solution  for  c  is 

c  =  — (l  -  cos [/gu(t-T) ] ) 

Therefore,  the  expression  for  f  reduces  to 


e   sf  =\ 


ujy 


u  n 


u  n 


sinh(un)  +  ^——{1  -  cos[/gU(t-T)]} 
u  u  ^ 


u|y 


sinh(uy}  +  ^i_i{l  -  cos[/gu(t-T)]} 
u        u 


(y<n) 
(y>n) 
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Inverting  to  obtain  G  (x  ,y ;  £  ,r| ;  t ,  t) 


TTjo         U 


sinh  (uri)  cos  [u  (x-£)  ]  du         (y    <    n) 


G    =< 


e 


+ 


u 


sinh (uy) cos [u (x-5) ] du         (y    >    n ) 


,«>      u(y+n) 
+   —        [1    -    cos/gu(t-T)]    cos    u(x-Odu 

7T  I  o  U 


Recalling  the  known  integral 


00   -ax   -a*x 
e   -e 

n       X 


cos (mx) dx=log 


/a'2  +  m2 
/a2  +  m2 


then  G  may  be  expressed  as 


G(x,y;£,n;t,T)  = 


1  , 
t—  log 


/(x-£)2  +  (y-n) 2 

.  /(x-?)2  +  (y+n)2 


-1-  eu(y+n)     U    ~    cos/^I(t-T)3    cos[u(x.g]du 

7T  J  0  U 


G  (x,y ;  t,  ,r\ ;  t,  t)  -  Fluid  of  Finite  Depth  h. 
G  must  satisfy  the  following  governing  equations: 
Zn   _  32G  ^  32G  _  „„  „*/„_„*  *_  -h  <  y  <  0 


VZG 


3X2  +  ^T  =  <5(x-Q6(y-ri)  for 


-<»   <   X   <   +°° 


G.  ,  +  gG   =0    on  y  =  0 ,  the  free  surface 

^  boundary  condition 


G   =  0 

y 


G  =  Gt  =  0 


on  y  =  -h,  the  bottom  boundary 
condition 

for  t=T,y=0,  the  initial 

conditions 


59 


Substituting  for  G  its  Fourier  transform  representation  f, 
the  governing  equations  become 


1)  eiu?(fyy  -  u2f)  =  6(y  -  n) 


2)  f . .  +  gf   =0    on  y  =  0 
tt    3  y  J 


3)  fy  =  0 

4)  f  =  f   =  0 


on  y  =  -h 
for  t  =  t ,  y 


=  0 


For  y  <  n /  the  solution  to  eq.  1)  subject  to  eq.  3)  is 


e        f  =   a  cosh[u(y+h)] 


For  y  >  n /  the  solution  to  eq.  1)  subject  to  eq.  2)  is 


e    f  =  c  cosh  (uy)  -  sinh(uy) 


Require  f  to  be  continuous  at  y  =  n  and  that  the 
derivative  of  f  have  a  discontinuity  of  magnitude  -1, 
that  is 

y=n+       J Jy=n- 


=  -l 


From  the  continuity  requirement  above 


1UU 

e   ^f 


cosh{u (y+h) } [c  cosh(un)  - 
cosh{u (n+h) } [c  cosh(uy)  - 


_tt 

ctt 
gu 


sinh (un) ] 
sinh (uy) ] 


(y  <  n) 
(y  >  n) 


From  the  requirement  on  the  derivative  of  f 


c^   +  gu  tanh(uh)c  =  ,  ,  ,  > 

tt    3  cosh(uh) 


Subject  to  c  =  c   =  0    for  t  =  x 


The  solution  for  c  is 
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u  sinh(uh) 
where  a  =  /gu  tanh(uh) 


{1  -  cos [a (t-x) ] } 


Therefore,  the  expression  for  f  reduces  to 


eiu?f 


sinh (un)   . 
u  cosh(uh) 


cosh  [u  (y-i-h)  ] 

cosh[u  (n+h)  ]  {1  -  cos [a (t-x)]} 


(y  <  n) 


u  sinh (uh) cosh (uh) 

i  r  /  ,  ■  i_  \  i      sinh  (uy) 

cosh  u  (n+h)  J   -  '  t*!v    + 

|_  u  cosh(uh) 

cosh [u (y+h) ] {1  -  cos[o(t-i)  ]  }" 


u  sinh (uh) cosh (uh) 


(y  >  n) 


Inverting  to  obtain  G  (x,y  ;  E,  ,r\  ;  t ,  x) 

cosh [u (y+h) ] sinh (un)     r  ,   >-\-ij   ,   „   \ 
l—  *eL,„u\ — ■ — —  cos[u(x-£)  ]du  (y  <  n) 


o 


u   cosh(uh) 


G    =    { 


cosh   uri+h      sinh   uy  r     ,      rN1.-,       /      .    _\ 

—    /  7  .  x —  *   cos  [u  (x-£j  ]du  (y  >  n) 

0      u  cosh (uh)  ^      2 


cosh  [u  (n+h)  ]  cosh  [u  (y+h)  ]  [l-cos{a  (t-x)  }]__,_  r„  ,__  rx  1/q„ 

_____ —-j a _ . — . — _— cos  iu^x—c;  JU.U 

.        coshMuh)  u  tanh  (uh) 
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Recalling  the  useful  identity- 
cosh  [u (y+h) ]  =  e  -^cosh(uh)  -  e    sinh(uy) 
then  G  can  be  expressed  in  the  following  form 

1  f   e  ^sinh(uri) 


u 


-cos  [u  (x-£)  ]  du      (y  <  n) 


G  =  < 


e   smh(uy) 
u 


cos [u (x-£) ]du     (y  >  n) 


■°°   -uh 


1  f   e    sinh  uy  sinh  uri      ,  ,   rxn^ 

+  —     t  i    u  \ — * — ~  cos  [u  (x-^)  ]du  + 

it  J  0      u  cosh  (uh)  s 


+  ± 

IT 


•rocosh[u(y+h)3cosh[u(n+h)]>[l-cos{a(t-T)}]CQs[u(x_U]du 
n       cosh* (uh)  u  tanh(uh;  *' 


where  a  =  /gu  tanh(uh) 

Using  the  known  integral  expression  for  the  log  given 
previously,  G  can  finally  be  expressed  as 


G(x,y ;£, n; t,i) 


~—  log 


/(x-c)2  +  (y-n)2 

/(x-£)2  +  (y+n)2 


+  i 

7T 


00   -uh  . 


J  o 


e    sinh  (uri)  sinh  (uy) 
u  cosh (uh) 


cos [u (x-£) ] du  + 


+  i 


'°°cosh  [u  (n+h)  ]  cosh  [u  (y+h)  3     [l-cos{a  (t-i)  }] .r„,„_n  ,  A„ 

cosh2(uh)  u   tanh(uh)         coslulx    UJ 


J  0 


where   0    =    /gu   tanh(uh) 
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APPENDIX  3 
THE  UNSTEADY  WAVEMAKER  SOLUTION 

Consider  a  semi-infinite  body  of  incompressible, 
inviscid,  fluid  having  a  free  surface  and  infinite  depth. 
Let  the  limiting  boundary  execute  infinitesimal  displace- 
ments from  an  initial  vertical  plane  position  given  by 
s  =  SF(y,t)  where  S  has  dimensions  of  length  and  therefore 
F(y,t)  is  dimensionless .   Initially  F(y,0+)  =  0.   It  is 
assumed  that  F(y,t)  is  sufficiently  smooth  and  vanishes 
sufficiently  rapidly  as  y  approaches  -°°.   The  problem  is 
described  mathematically  in  terms  of  a  velocity  potential 
as  follows : 

V2*  =  0  +  0=  0    y  <  0,  x  >  0 

subject  to   cj)    +  g<j)   =0   on  y  =  0 ,  the  free  surface 

y  boundary  condition 

4>  y    0  as  y  >    -°° 

cj>x  =  SF(y,t)    on  x  =  0 

4>   =  <J>t  =  0     for  t  =  0+;  y  =  0,  the  initial 
y  conditions 

The  relevant  time-dependent  Green's  Function  asso- 
ciated with  the  above  problem  is  easily  obtained  by  the 
method  of  images  and  is  given  by: 
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G(x,y;£,n;t,T)    =    -   ^7  lo<3 


/(x-g)  *+(y-n)  */Tx+g)  z+(y-n)  * 
_/(x-S)  2+(y+n)  2/(x+£)  2+(y+n)  2  _ 


+ 


+  i  eu(^+n)     [1    -   cos{/^T(t-T)}]     [cos{u(x_a}    f 

7T      J0  U 


+  cos{u(x+£) }]du 


Note:   The  only  singularity  in  the  region  occupied  by  the 
fluid  is  at  (x,y)  =  {E,  ,r\)  .      All  other  singular 
points  are  outside  the  boundary  enclosing  the  fluid 

Applying  Green's  Theorem  of  vector  analysis  (in 
2-dimensional  form)  to  the  region  occupied  by  the  fluid 
exclusive  of  the  singular  point  (£,ri)  one  obtains: 

Green's  Theorem 


y 


-►  x 


R 


(<})1V24)2  -  *2V24>1) 


dA  = 


R    Urn) 

0 


8^2  3^t 

(J)  (4), -5--  -  c{)„-7— )ds 

1  dn  2  dn 
c 


Let   (J>1  =  G(x,y;C,n;t,T) 
<j>2  =  <j>t(x,y;t) 

Then  substituting  into  the  above  expression  for  Green's 
Theorem  , 


(GV2c})t  -  4>tV2G)  dA   =    0    (G(J>        -  <j>  ,G    )  ds    +    0       (G<J>t     -  *tGn>  ds 
R'  •'c  '  ce 

where    R'    =    R  less    the   singular   point    (£,ri);    ce    is    a   circle  of 

vanishingly    small    radius    enclosing    {€,  ,r\)  . 
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'2  a     - 


However,    V*<J>   =    0    implies    V    <$>,    =    0    and    V   G   =    0    except   at 
(£  ,n)  .      Therefore 

(G<j>.       -    <j>.G  Jds    =   0       (G<j>.       -    <j>.G    )ds 


tn 


t  n 


tn 


t  n' 


On  ce/  ds  =  rd6,  3/8n  =  -  j-   where  r  =  /(x-£) 2  +  (y  -n) 2 


With  this  substitution  the  above  expression  reduces  to 

•2tt 


d>    (Gcj>.       -    (j>.G    )    ds    =    ll]JJ 
rtn        rt   n  r+0 

c 


-(G<j>.   -  4>G  )  rd9 
rtr    r  r 


cf>   is  regular  throughout  R;  therefore,  the  only  contribu- 
tions to  the  above  limit  must  come  from  G  and  G   (assuming 

r         3 

the  limit  of  the  integral  is  equal  to  the  integral  of  the 
limit)  and  are  given  by 


, ,  x  lim  r    _,,  „  lim 
(1)  r-*0  [rGl  "  r>0 


log  r 

2tt 


r  ►()     r     r  *•() 


2iir 


2n 


Therefore 


*  (G*tn  "  Wds 
c 


-(f).  U,n;t)  f2TT 


2tt 


de  = 


-*tU,n;t) 


Thus 


<J>t(C,n;t)    =  <j>    (4>tGn  -   G<j>tn)ds 
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Integrating  the  above  expression  from  0  to  T 


<j> .  (§,n;t)dt  = 

o    t 


t    r 


o     L 


(<J>.G      -    G<J>.     )ds 
rt   n  Ttn 


dt 


Note  that  from  the  behavior  of  G  and  G   at  °»  and  from  the 

n 

• 

prescribed  boundary  conditions  it  is  easily  shown  that  0 

J  c 
reduces    to 

o 


0 


+ X  +  + 

y=-oo       ^J  x=+°°  ^  y=0 


x=0 


Onx=0,  9/3n  =  -3/3x.   On  y=0,  9/3n  =  +3/9y.   Therefore 

]  <J>t(C/n;x)  =  <KS/n;T)  -  4»<4#n;0)  = 


((KG      -    G4.     )dx    -  (<|>.G      -    G<J>.     )dy 

Jo         Jot   y  Yty  jo         t   x  rtx 


dt   = 


(y-0) 


(x  =  0) 


(G<j>.        -    <J>.G    )dt 
J      v    Tty         yt    y' 


(y-0)  L 


.—00 

t 

dx    + 

J 
(x- 

0 

=  0) 

• 

fX 


<G*tx   -   *tGx>dt 


dy 


'0 

(y  =  0) 


(GcJ>      +   -V1) 

y  g 


fX 


'tt. 


G.  (<J>      +   -=^)dt 

o    t   Ty  g 


dx    + 


+ 


Jo 
(x  =  0) 


(G«J)x    -    <|.Gx) 


(GJ       -    ()iGt     )dt 
v    t    x  tx 


dy 


For  the  problem  at  hand 


G  (0,y;£,n;t,T)  =  G.v(0,y;5,n;t,T)  =  0 


x 


<f>x(0,y;0)  =  0 


tx 


Ty   g 


tt 


0   on  y  =  0 
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G(x,0;£,n;T,T)  =  Gt(x,0;C,ri;T,T)  =  0 

<j>  (x,0;0)  =  <j>t(x,0;0)  =  0 

(J)(x/y;0)  =  0   (as  the  fluid  is  assumed  initially 

at  rest) 


Therefore 


c|>  ( £  /  n ;  t  )  = 


[G(0,y;£,n;T,T)<b  (0,y;x) 
J  o  x 

(x-0) 


-  rGt*xdt]dy 


Let 


where 


s  = 


v 


S  eUoy  sin(wt)  H(t) 

wS  eU°y  cos(wt)  H(t)  =  <j)  (0,y;t) 


x 


=  -w2S  eUoy  sin(wt)  H(t)  + 

+  wS  eU°y  cos(wt)6(t)  =  «|>   (0,y;t) 


u   =  w  /g 
o      /  3 

H(t)  is  the  Ileaviside  step  function 
6 (t)  is  the  Dirac  Delta  function 
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Substituting  these  expressions  into  <J>(£/n;f)  above,  it  is 
easily  shown  that 


<J>  ( £  /  n ;  t  )  = 


-wS  cos(wt)H(t) 

2tt 


%y  ma  I  liliYznl 


e  °*  log  |  fTTTyTnt^  /  dy  + 


-2wS     e   cos(ug) 


lo      u 


L  J 


(u+u  )y  , 

e     o  ■*  dy 


rT 


'gu  sin  [/gu(t-x)  ]  cos  (wt)H  (t)  dt 


i-  a 


du 


Of  more  direct  interest,  however,  are  the  waves 
in  lieu  of  the  velocity  potential.   The  expression  for 
the  waveheight  is  obtained  from  the  linearized  free  surface 
boundary  condition  given  by 


n(£/T)    = 


|   <J>tU,0;t) 


Thus 


C°S|U1V  cos{/gU(t-T)  }cos(wt)H(t)dt 

o       u   "*"    uo     L  Jo 


du 


or 


nU,i)  = 


Wb 


cos (uC) 


71  Jo     u  +  u0     L 


sin (wt) -sin (/gui)       sin  (wt) +sin ( /gur) 
w  -  /gu  w+  /gu 


du 


where   u^   =  w2/g 
o  '  3 
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APPENDIX  C 
ASYMPTOTIC  BEHAVIOR  OF  CERTAIN  INTEGRALS 


1.    The  asymptotic  behavior  of  the  integral  given  below 


is  desired. 


C  d*-l 


a2    i(a2x-at)  , 
-  e         da 


The  contour  C-,  has  been  chosen  so  that  it  corresponds  to 
the  steepest  descent  path  from  the  origin  to  infinity. 
Introduce  the  change  in  variable  ii  =  a2 (x/t)  -a.   With 
this  substitution,  the  complex  variable  t  is  real  along 
this  contour.   The  resulting  integral  expression  is  the 
following : 

•  00 

(o2/ok-l)    da/dx  e~tT  dx 

J  o 

The  asymptotic  behavior  of  this  integral  can  be  ascertained 
by  expanding  the  integrand  about  t  =  0  and  integrating  the 
resulting  expression  term  by  term.   However,  the  integrand 
is  regular  about  this  point  and  therefore  the  Taylor  series 
expansion  is  applicable.   The  series  expansion  is  given 
below  and  is  readily  obtained  through  repeated  differentia- 
tion of  the  integrand. 


69 


oz      da 


"♦_!   dT 


a"-l 


=  f(T)   =  f(0)  +  f'(0)x  +  f"(0)T  + 


=  0  +  0  -  ix2  + 


The  first  nonzero  term  of  the  Taylor  expansion  is  of 
order  0(t2).   The  dominant  term  in  the  asymptotic  expan- 


sion of 


is  then  proportional  to: 


o2    i(a2x-ot)  , 

e  da  =  -l 


a»-l 


2   --tT 


Jo 


x   e    dx  +  . . . 


Introducing  the  change  in  variable  tx  =  y,  the  above  ex- 
pression can  be  expressed  as  follows: 


a2    i(a2x-at)  ..      >x.-3 
e         da  =  -it 


a--l 


y2e  ^dy  +  . . . 


=  -it~ir(3)  +  .  .  . 


The  dominant  term  in  the  asymptotic  expansion  of  this 
integral  is  therefore  of  order  0(t~3). 


Consider  the  following  integrals 

i (a2x+at) 


Ji  ■ 


h  = 


X4    " 


f(o)    e" 


dcr 


f(olei(o2x+t)    da 

o        a 

°°  1M    ei(a2x-t)    da 


where      f(a)    =   a   /a    -1 
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The  asymptotic  behavior  of  these  integrals  can  be  obtained 
by  deforming  the  path  of  integration  in  the  complex  a  plane 

along  a  contour  with  the  general  features  of  C,  below. 

icf.^ 

cr-plane 


A  particularly  convenient  equivalent  contour  for  the  evalu- 
ation of  I,,  I3,  and  I.    consists  of  the  steepest  descents 
path  from  the  origin  to  infinity  and  then  along  the  arc 
Re   to  the  real  axis.   Such  a  contour  is  shown  below. 

where  C^  is  the  steepest 
descents  path 

C  =  Re19,  R  ■*  °° 


For  integrals  I.,  and  I.,  the  steepest  descent  path  is 

iTr/4 
simply  the  ray  a   =   r  e     .   For  integral  I, ,  the  steepest 

descent  path  is  not  as  easily  described  but  no  further 

complication  is  added  to  the  analysis. 

As  the  integrand  is  in  all  cases  regular  at  a  =  0 ,  it  may 
be  expanded  in  a  Taylor  series  expansion  about  a  =  0  and 
the  resulting  expression  integrated  term  by  term.  By  in- 
spection, the  first  two  terms  of  the  Taylor  series  are 
equal  to  zero  for  I, .   Only  the  first  term  is  zero  for 
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integrals  I~  and  I..   The  contribution  to  the  integral  is 
then  of  the  order  0(t~3)  for  integral  I-,  and  of  order 
0(t~")  for  integrals  I~  and  I.. 


1 


The  asymptotic  behavior  of  integrals  I,,  1~,    and  I.  can  be 
described  as  below: 


o 
I,  ^  i-rr  Res(l)  +  i/*  +    0(t-3) 

(R->-°°) 


(^)ei(x+t)  +  0(t"3) 


I3  ^  i-n  Res  (1)  + 


+  0(f)  = 


(ijL)ei(x+t)  +  0(t-2) 


R 
(R->-°°) 


0 

I4  ^  iTT  Res(l)  +  y^    +    0(t~2) 

(R-+oo) 


(^L)ei(x-t')  +  0(t"2) 
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APPENDIX  D 
LAURENT  EXPANSIONS 


To  obtain  an  asymptotic  expansion  for     ,  it  is  con- 


venient  to  obtain  the  Laurent  expansion  of 


o2   da 


a"-l  dx 

about  x  =  0  such  that  the  effect  of  the  pole  at  a  =  1  has 
been  subtracted  off  with  the  remainder  of  the  expression 
regular  at  a  =  1 •   Algebraic  difficulties  are  minimized  if 
a2/a'+-l  is  first  expanded  by  the  method  of  partial  frac- 
tions.  It  is  easily  shown  that 


aH-l    4 


A-  +  r±.  +  =!_  +  _L 


(t+i    o-i    a+1        a-1 


do 


Consider  first,  the  expansion  of  [o''/o''-l]  — = — 

1  +   +     dx 


o;  do  + 


o^-l  dx 


-l 


rr  +  rr-  + 

l+2iv)  +  yx  /2    (l-2iv)  +  yt 


+ xj-  +    rj- 

(l  +  2v)  +  yx  /2    (l-2v)  +  yx /2 


Y 


8x 


Tz 


This  expression  is  easily  placed  in  the  form 


2      da 


+  +         -y/8 


a?-l    di  l-2v 


-y2 


t/2    +    (iz2v; 


-1 


l-2v         l+2iv  V2 

+     T 
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-i 


-1 


i-2iv   +    T  %         l+2v   +       y2  _ 


=  -y/8 

l-2v 


Ty2  +  (iZ2v) 

Y 


-1 


+  x 


-y2 


l-2v 


{l-4v   -    4v2}t     /2         y{ (l+4v2) 2    -    8v(l+2v) 2) 


(l+2v)(l  +  4v2) 


(l+2v) A  (l+4v^) 


2  \   2 


+    -Y2((l  +  4v2)3    -    4v(l+2v)  3(3-4v2)  }t  /2 
(l+2v)  3     (l+4v2)  3 


+    .  .  . 


da 


2.  The  expansion  of  [o_/o_-l]—^ —  is  obtained  in  a  similar 


manner. 


a2   da 


o4-l  di 


-l 


T7-    +     = T7-    + 

(l  +  2iv)  -  yt/2    (l-2iv)  -  yt  /2 


-1 


—17-  +  

(l+2v)  -  yt      (l-2v)  -  yr 


Yz 


8t 


-Y 


3y    inspection,  the  expansion  of  this  function  is  identical 
to  tne  one  previously  obtained  except  for  sign  changes  on 
alternate  terms.   The  resulting  expression  is  readily 


placed  in  the  following  form: 
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a2   da 


v8 


a"-l  dx    l-2v  *- 


y2  _  (i-2V) 

Y   J 


-|-1 


+  -J- 


Y 


-t  n    +  {1-4v-4v2}t"/2    y{(l+4v2)  2  -  8v(l+2v)  2} 


8  L  i_2v     (l+2v)  (l+4v2) 


(l+2v)2   (l+4v2) 


2  \  2 


y2{(l+4v2)3  -  4v(l+2v)  3(3-4v2)  H^2 

•    ■  -  -  —  .  —  .—.  —  .ii,  -^-   •  •  • 

(l+2v) 3  (l+4v2)3 


3.  For  values  of  v  =  x/t  <<  1/2  or  v  =  x/t  >> 1/2  there  is  no 
particular  advantage  gained  in  subtracting  off  the  effect 
of  the  pole  at  o   =   1.   As  before,  algebraic  difficulties 
are  minimized  by  expressing  a2/ a1*-!   as  a  sum  of  linear 
terms . 


a2      da 

i      +   "i      +   "I      +      1 

1  da 

ok-l   dx 

a+i        a-i        a+1        a-1 

4   dx 

Using  the  familiar  expansion 


=  1  +  x  +  x2  +  x3  +  x1*  +  x5  + 


1+x 


the  following  expressions  may  be  obtained: 


f±  !!±  =    2V3/2    e^/4 


a^-1     dx 


-Vz  4V1/2i1/2(l+16V,+  ) 


1-16V' 


(1-16V14) 


4  \   2 


2\   2  ,    //it72  s    k  \  ,,  /l 


+    4Vi{3  +  12(4V2)2-K4V2)'t}T 


«♦  \   3 


(1-16V) 
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a2      da 


,2v3/2    eiTT/4 


a1*-!    dx 


y2 


y2  4  y2 


+    4V/2i  /2  (1+16VH)    + 
1-16V1*  (1-16V1*)2 


,.  /2 


4Vi{3+12(4V2) 2+(4V2)"}x 
+    +    0  (  t) 


(1-16V) 


•♦  N    3 


4.         For   v   =    1/2,    the   following   Laurent   expansions    apply 


+  +  2 

=    +    Y 

a  1-1   dx 


Y"2t_1  1       -i    -Vz  3 

J +    v        X  -     + 


8  16 


32 


3    „J/i    .^1^.1^       .3T3A    + 


+    64    ^T 


128    T         256 


a2     da 


a4-l    dx 


]Cici+iY-iT-v.  +  ^  + 


V 


8  16 


32 


♦     J-VT1/*     +-JClx     --£-    Y3T/2     + 


64 


128  25b 
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APPENDIX  E 

ALTERNATIVE  EXPANSIONS  FOR  n(x,t) 


For  the  cases  where  V  =  x/t  <<l/2,  V  >>  1/2,  and 

V  =  1/2,  alternative  expressions  for  the  waveheight  may  be 
developed.   The  effect  of  the  simple  pole  at  a  =  1  is 
either  negligible  or  an  easily  evaluated  residue  contribu- 
tion in  these  cases.   These  alternative  expansions  for 
H(x/t)  have  intrinsic  merit  and  in  addition  provide  a 
check  on  the  uniformly  valid  expression  (35) .   Equations 
(6),  (10)  and  (24)  may  be  combined  and  rewritten  as  follows 

n(x,t)  =  Im  fe  ei(x_t)H(l-2x/t)  -  j    +  0(t"3)|   (vyi/2) 

C2 

The  asymptotic  behavior  of  n(x,t)  can  be  evaluated  for 

V  <<  1/2  and  V  >>  1/2  by  investigating  the  behavior  of 

/„  .   This  integral  (after  a  suitable  change  in  variable) 
c2 

can  be  expressed  as  in  (12) .   Laurent  expansions  of  the 
integrand  of  (12)  subject  to  the  restrictions  on  V  above 
have  been  developed  in  section  3  of  Appendix  D.   Substitu- 
tion of  these  expansions  into  (12)  yields  the  following 
result : 


2   it(a*V-a)do  *  4vVi(^  ~  ^  [~  -±-^—    + 


a 
e 


C2a*-1 


4iV{3+12(4V2)2+(4V2)'t}   V2    nt'/is 
(1  -  16V  )  3  ^ 


c00 

0 

r,  * 

_(1-16V'+) 

l)~ 

e        dx 

M 
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Integrating  term  by  term,  the  following  asymptotic 
expansion  for  /    is  obtained: 


3/ 

%  4V/2  e 


i(4  -  x)  T  r(i)t_1/ 


4V   4 


(1-16V*) 


4iV{3+12(4V2)2-K4V2)-}r(3)t-3/2  +  0(t-5/2) 

(1-16V")3       2 


(V>>l/2,V«l/2) 


The   waveheight    n(x,t)    can   now  be   written: 


Im    -j    ij-  ei(x    fc)    H(l-2x/t)    + 


■   /    t  TTv 

3A       "1(4V    "    4^ 
+     (-4)V/2    e         4V         4 


r(Lt"1/2  3/ 

— - +   0(t    /2) 

_  (1-16V4) 


(V^l/2) 


For  the  case  V  =  x/t  =  1/2,  the  Laurent  expansions 
given  in  section  4,  Appendix  D,  apply.   The  simple  pole  at 
o  =  1  is  now  coincident  with  tiie  saddle  point  at  a  =  1/2V. 
Therefore,  the  contour  C„  must  be  indented  and  the  residue 
contribution  accounted  for.   The  waveheight  n(x,t)  is 
given  by: 


itt   l  (x-t) 
n(x,t)=I  <    — r-  e 
m    4 


+  0(t  3) 


(V=l/2) 
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Substitution  of  the  appropriate  Laurent  expansions 
into  the  integrand  of  Jr      yields: 


^   e 


1(2 


-) 


72    -y2         3S2±T%    +    9/2    3/2  (    5/2" 

8  T  16  32  ulT      ; 


e         dx 


%   e 


1(2   "   T> 


"£r  <§)  t"  *  +  3jgir  ,3,  t-  3A  +  ^r  <§)  t"  *«  +  o  < t" ?/2 )' 


The  waveheight  n(x,t)  can  now  be  expressed  as 


iiT   i(x-t) 
n(x,t)  *  lm   ^  e       + 


+  (-)e 


"1(2  "  4> 


~  r(j)  f1/2  +  o(t"3/2) 


(V  =  1/2) 
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APPENDIX  F 
EXPERIMENTAL  DATA 

During  the  experiments  continuous  readings  of  wave 
height  and  time  were  made  at  two  locations  in  the  Ship  Model 
Towing  Tank  with  the  use  of  a  Sanborn,  Series  7700,  strip 
chart  recorder.'  The  motions  of  the  wavemaker  paddle  were 
also  continuously  monitored  and  recorded  simultaneously  with 
the  wave  data.   Time  was  measured  from  the  instant  of 
initial  motion  of  the  wavemaker  paddle. 

Reproductions  of  the  wave  records  obtained  as  the 
output  of  the  strip  chart  recorder  are  presented  in  order  of 
increasing  wavemaker  frequency.   Time  is  measured  horizon- 
tally and  wave  height  vertically  according  to  the  scales 
indicated. 
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APPENDIX  G 
COMPUTER  PROGRAMS 

The  facilities  of  the  Computation  Center  of  the 
Massachusetts  Institute  of  Technology  were  used  extensively 
in  the  comparison  of  theoretical  to  experimental  results. 
The  following  programs  were  developed  to  calculate  and  plot 
the  waves  predicted  by  equation  (35)  for  r\  (x,t)  : 

1)  Timewise  distribution  of  wave  height  for 
fixed  distance. 

2)  Spacewise  distribution  of  wave  height  for 
fixed  time. 

3)  Function  subprogram  for  the  error  function 
complement  for  complex  arguments. 
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APPENDIX  H 
NUMERICAL  ANALYSIS  OF  n(x,t) 

1 .  Short  Time  Solution. 

The  asymptotic  behavior  of  equation  (5)  for  the  wave- 
height  n(x,t)  is  readily  evaluated  for  large  values  of  t 
(or  x) .   This  behavior  is  described  by  equation  (35) .   There 
appears  to  be  no  such  expression  for  small  values  of  t 
(or  x) .   For  this  case  numerical  methods  must  be  used  with 
the  aid  of  computational  facilities  a  practical  necessity. 
For  very  small  values  of  t,  equation  (2)  may  be  evaluated  in 
its  present  form  -  integrating  numerically;  e.g.,  using 
Simpson's  rule.   As  t  (or  x)  grows  in  magnitude,  however, 
this  approach  becomes  unsatisfactory  due  to  the  increasingly 
oscillatory  nature  of  the  integrand  of  (2) . 

2 .  Programming  Difficulties  and  the  Solution. 

The  numerical  integration  difficulties  posed  by  the 
highly  oscillatory  nature  of  tne  integrand  of  (2)  can  be 
avoided  by  employing  the  contour  integration  methods  and 
appropriate  changes  of  variable  used  in  determing  the 
asymptotic  oehavior  of  equation  (5).   Here  tne  integrals  I,, 
I~,  I,  and  I.     (defined  by  (5))  are  regarded  as  contour  in- 
tegrals in  the  complex  o  plane  and  equivalent  contours 
which  facilitate  their  evaluation  have  been  determined. 
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Integrals  I,  ,  I .,  and  I4  can  therefore  be  represented 
by  (see  Appendix  C) : 


I,  =  iTTRes(l)  + 


CSD  «*"1 


a2    it(a2V+a), 
e         da 


I3  =  iTTRes(l)  +  e 


it 


a    1a  x  , 
e     da 


CSD  ""-1 


I   =  iiTRes  (1)  +  e 

4 


-it 


csD  "*-1 


a    la  x  . 
e    da 


Introducing  the  change  of  variable  ix   =  o   V+o    into  I, 

i  i 

and  ix   =  a2  into  I->  and  lA,    these  integrals  can  be  written 

2  -5  4t 


I,  =  iiTRes(l)  + 


a2   da    -tx,  , 


ax 


o  aH-l  dx 


I3  =  iirRes(l)  +  e 


it 


,00 

a   da    -xx 


e  '  2  dx 


o  a  -1  dx 


I4  =  iTTRes(l)  +  e 


.00 

-it     a    da 


-xx, 


e    2dx 
o  a1*-!  dx  2 


When  expressed  in  these  forms,  the  oscillatory  nature  of 
the  integrand  has  been  removed  and  in  fact  it  now  becomes 


exponentially  small  as  x   and  x   become  large.   Numerical 

1  2 
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evaluation  of  these  integrals  can  now  be  accomplished  by 
a  variety  of  methods  for  small  and  large  t  (or  x) . 

A  similar  method  can  be  used  to  evaluate  I~.   Its 
representation  is  given  by  equations  (9)  and  (23)  for  the 
cases  V  =  x/t  <  1/2  and  V  >  1/2.   Considering  only  equation 
(23) ,  I~  can  be  written: 


/-is   f    o2        it(a2V-a),    f 
I0  =  17T Res  (1)  +  I    e         da  + 

2  Jr   n--l  h 


■'C1  ak-l  JC2  ok-l 


eit(a2V-a)da 


(V  >  1/2) 


where  C,  is  the  steepest  descent  path  from  the  origin  to 

infinity  and  C-  includes  both  steepest  descent  paths  from 

the  saddle  point  to  infinity.   Introducing  the  change  in 

variable  ix   =  o2V-o    into  /    and  -1/4V  +  ix,  =  o2V-o    into 

j       ,    I„  can  be  expressed  as: 
c2    z 


I2   =    iTTRes(l)    + 


-it/4V 
+   e        ' 


1  \J>J  r\ 

o'      da         -ti,      , 
e        l    di2    + 

o    a1*-!   dij 


af__  da_  e   tx2    dT^       (v   >    1/2) 


o    a4-l   dx2 


A  similar  result  can  be  obtained  from  equation  (9)  and  the 
case  V  <  1/2  -  the  only  difference  being  a  change  in  sign  of 
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the  residue  term.   Numerical  evaluation  of  I~  is  now 
readily  accomplished.   Finally,  equation  (5)  which  relates 
integrals  I,  ,  I~,  I-.  and  I.  to  the  waveheight,  is  used  to 
complete  the  numerical  evaluation  of  ri(x,t). 
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